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ABSTRACT 

The active galaxy Markarian 421 underwent a substantial outburst in early 
2001. Between January and May of that year, the STACEE detector was used 
to observe the source in 7-rays between the energies of 50 and 500 GeV. These 
observations represent the lowest energy 7-ray detection of this outburst by a 
ground-based experiment. Here we present results from these observations, which 
indicate an average integral 7-ray flux of (8.0 ±0.7± 1.5) x 10~ 10 cm _2 s _1 above 
140 GeV. We also present a light curve for Markarian 421 as observed by STACEE 
from March to May, and compare our temporal, as well as spectral, measurements 
to those of other experiments. 

Subject headings: galaxies: active — BL Lacertae objects: individual (Markarian 
421) — gamma rays: observations 
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1. Introduction 

The blazar Markarian 421 is one of only a handful of astrophysical objects detected in 
the very high energy (VHE) 7-ray regime, between 10 GeV and 100 TeV. At a red-shift of 
0.031, it is the closest BL Lac object seen by EGRET, and the first to be detected by a 
ground-based TeV instrument (Punch et al. 1992; Petry et al. 1996). Along with Markarian 
501 (z=0.034), it is one of the most prominent sources of extra-galactic VHE 7-rays, and has 
been regularly monitored by atmospheric Cherenkov telescopes since its TeV identification 
(see Catanese & Weekes 1999, for a review). 

Markarian 421 is an X-ray selected BL Lac object, and as such, it exhibits the double- 
humped spectrum that is characteristic of blazars (e.g. Takahashi et al. 2000). The lower 
energy hump, believed to arise from synchrotron radiation, peaks in the X-ray range, while 
the higher energy hump, generally attributed to inverse Compton (IC) scattering of soft 
photons, peaks in the GeV to TeV energy range. 

The details of the emission processes that produce the IC hump are, as yet, unresolved. 
For example, the origin of the soft photons that seed the IC component of the spectrum 
is an outstanding question in the study of blazars. Although synchrotron self-Compton 
models, in which the seed photons originate from synchrotron radiation within the blazar 
jets, are generally favored for Markarian 421 (Coppi 1992; Coppi & Aharonian 1999), other 
competing models exist. These include external Compton models (Dermer, Schlickeiser, & 
Mastichiadis 1992; Sikora, Begelman, & Rees 1994), proton-induced cascades (Mannheim 
1993), and proton synchrotron models (Aharonian 2000; Miicke & Protheroe 2001). 

Markarian 421 has been given to strong periods of flaring activity over the past few 
years. Of note was a flare in 1996 during which the source was detected at more than 
ten times the flux of the Crab Nebula (Gaidos et al. 1996). In the early part of 2001, a 
rather impressive flare was again observed (Boerst et al. 2001), with reported TeV fluxes of 
comparable magnitude to those recorded in 1996, but of much longer duration (Krennrich 
et al. 2001). 

It was during the 2001 flare that STACEE-48, an intermediate incarnation of STACEE 
(the Solar Tower Atmospheric Cherenkov Effect Experiment), was commencing operation. 
STACEE was able to observe the activities of Markarian 421 for much of the flaring period. 
These observations are of particular note as they represent the only 7-ray detection below 
200 GeV during this flare. 
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2. The STACEE-48 Detector 

STACEE is a low threshold ground-based 7-ray detector located at the National Solar 
Thermal Test Facility near Albuquerque, New Mexico. STACEE uses the large steerable 
mirrors (heliostats) of an existing solar research facility to collect Cherenkov photons from 
particle cascades in the atmosphere. The photons are focused, via secondary optics lo- 
cated on a central tower, onto an array of photomultiplier tubes, and their arrival times 
are recorded and processed by high-speed electronics. Because of their very large collection 
areas, wavefront sampling detectors, such as STACEE and CELESTE (de Naurois et al. 
2002), are currently the only ground-based instruments capable of detecting 7-rays below 
200 GeV. STACEE was commissioned for its first science observations with 32 heliostats 
(STACEE-32, see Hanna et al. 2002) and was used to detect 7-ray emission from the Crab 
Nebula in 1998 and 1999 (Oser et al. 2001). 

STACEE-48 (Covault et al. 2001) was an upgrade of the prior 32-channel experiment in 
which 16 new heliostats were added, bringing the total to 48, each with a mirror area of 37 
m 2 . The introduction of these new heliostats had the dual advantage of increasing the total 
collection well as producing a more favorable detector geometry. In addition, by 

incorporating these new channels into the trigger, the energy threshold of the experiment was 
lowered from 190 GeV to 140 GeV due to improved suppression of the night sky background. 

The trigger hardware for STACEE-48 was re-instrumented with custom built modules 
that handle the tasks of signal delay and trigger formation (Martin & Ragan 2000). The 
delays are constructed from a 125 MHz pipeline capable of storing and delaying one hit per 
channel in each 8 ns bin. The system has a coincidence resolution of 1 ns, by virtue of a 
vernier encoding scheme, and a maximum possible delay of 2048 ns. It applies delays to 
each channel such that an idealized Cherenkov event will produce a coincident signal across 
the entire detector. Channels are grouped by geography into subclusters of eight heliostats 
each, and a trigger was formed if four of the six subclusters were hit within a 28 ns window. 
Individual subclusters registered a hit if five of their eight constituent channels fired within 
a 12 ns window. 

Improvements from STACEE-32 were also made on a number of other fronts. We 
aligned the mirror facets on all 48 heliostats using a laser look-back system. The phototubes 
were re-calibrated, and then re-distributed in the camera to better balance the detector 
response. We used the full Moon to align the optics from end to end, which resulted in 
important improvements to the optical throughput. We also reduced the background light 
reflected from the asphalt on the heliostat field by up to 40% by resurfacing the field with a 
dark sealant. 
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We have conducted extensive simulations of the STACEE-48 detector and its response 
to both 7-ray and hadronic events. We employ the CORSIKA air-shower simulation package 
(Heck et al. 1998) with the QGSJET hadronic interaction model (Kalmykov & Ostapchenko 
1993), together with custom optical and electronics simulation packages to obtain the detec- 
tor response. We have been successful in reproducing various detector characteristics with 
the simulations, including the detector rate at zenith and the trigger rate as a function of dis- 
criminator threshold. From simulations of our effective area convolved with a 7-ray spectrum 
of differential index a = 2.1 (see Fig. 1), the energy threshold for our observations (defined 
as the peak of the detected differential spectrum) was determined to be 140 ± 20 GeV. The 
error represents the shift in the peak location as the spectral index is varied between 1.7 and 
2.5. 



3. Observations 

STACEE-48 was used to observe Markarian 421 in the early part of 2001 (Boone 2002). 
STACEE uses an on-off technique for source observations as described in Oser et al. (2001). 
We take on-source data for 28 minutes, with a corresponding off-source run of 28 minutes to 
form an on-off pair. Our observations from March to May included about 26 nights, yielding 
78 on-off pairs from 36 hours of source observations and an equal amount of background 
observations. 

We have found that the subcluster firing rates, which are driven primarily by the night- 
sky background, are a powerful diagnostic of abnormal fluctuations in sky conditions. In 
good-quality pairs, the time evolution of the subcluster rates for both on and off-source 
runs are very linear, and well correlated. However, in runs with passing clouds or high 
haze, the rates are irregular and uncorrelated. Thus, removing runs whose subcluster rates 
exhibit deviations from a linear time-evolution has proven to be an effective quality selection 
criterion. Of the original 78 pairs, 17 were removed from the data set using this criterion. 

Three on-off pairs were removed because they showed abnormally high rates from back- 
ground light, and two pairs were removed due to heliostat malfunctions. Other pairs were 
trimmed to remove short-term contamination from an otherwise acceptable run. Remov- 
ing and trimming runs reduced the usable data to 56 on-off pairs, totaling 22.0 hours of 
observations. 

After applying all the selection criteria, the data set consists of 317, 045 ±563 on-source 
events and 307,641 ± 555 off-source events, yielding a 9,404 ± 790 event excess in the on- 
source data. This corresponds to a raw average excess of 7.1 ± 0.6 events/min, where the 
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quoted error is statistical only. 

Of particular concern to non-imaging atmospheric Cherenkov detectors like STACEE is 
the relatively bright star HD 95934 (magnitude 6.16 in the B band), which is within a few 
arc-minutes of Markarian 421, and well within the STACEE field of view. A calculation of 
the expected photon flux from this star and the trigger multiplicity required by the detector 
indicates that we should see much less than one false trigger per night from the star alone. 
However, random fluctuations in the excess noise generated on each channel by the starlight 
may promote an otherwise sub-threshold cosmic-ray event. Such promotions of sub-threshold 
background events would only occur in on-source observations, as there is no corresponding 
star in the off-source field of view. In this way, it may be possible to create a fake signal. 

To study this star effect, we spent some time observing HIP 80460 (RA=16.42 hr, 
Dec=37.39°), a star with a magnitude and declination comparable to HD 95934. We took 
12 on-off pairs, consisting of approximately five on-source hours. After applying analogous 
quality requirements to those described above, 11 pairs and 4.3 hours of data survive. These 
data indicate a raw excess of 1.9 ± 1.2 events/min from the star effect. 

After correcting for this fake signal, the detector deadtime, and the distributions of 
detector configurations (as described in the next section), we detect an average 7-ray rate 
from Markarian 421 of 7.7 ± 0.7 ±1.2 7/min, where the first error is the statistical error on 
the Markarian 421 measurement alone, and the second represents the statistical error in the 
observations of HIP 80460, used to correct for the effect of the star HD 95934. 

4. Flux Determination 

Although the source is expected to be highly variable, it is useful to calculate an av- 
erage flux over all the STACEE-48 observations. Detector simulations indicate that the 
trigger rate varies significantly as a function of zenith angle, and, to a lesser extent, detector 
thresholds. In order to calculate quantities that are relatively independent of these effects, 
we corrected the measured rate at a given detector zenith angle and threshold to the corre- 
sponding expected rate at a fiducial configuration (a zenith angle of 5°, near Markarian 421 
transit, and a detector threshold of ~ 5 photoelectrons on each phototube). To do this, we 
simulated the detector response to an assumed source spectrum at the fiducial configuration, 
as well as the observation configuration. We then scaled the measured rate by the ratio of 
these responses to obtain a corrected const ant- configuration curve. 

These corrections are sensitive to assumptions about the functional form of the incident 
7-ray flux. Thus, we must first choose a hypothesis for what the 7-ray spectrum might be in 
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our energy range. Figure 2 is a spectral energy distribution of Markarian 421 showing mea- 
surements by EGRET, averaged over five years of observations (Hartman et al. 1999), and 
Whipple during the 2001 flare (Krennrich et al. 2002). The STACEE butterfly was obtained 
by assuming a simple power-law spectrum with a given spectral index. The normalization 
was then fit such that the convolution of the assumed spectrum and the detector effective 
area matched the actual rate measured by STACEE. This procedure was done for a range 
of differential spectral indices from 2.00 to 2.20, yielding a derived differential photon flux of 
(6.4±0.5 stat ±1.3 SJ , s ) x 10~ 9 cm _2 s ^TeV" 1 at 140 GeV. Here the statistical error is the error 
on the Markarian 421 measurement. The systematic error is obtained by adding in quadra- 
ture the statistical error from the HIP 80460 measurement (1.0 x 10~ 9 ), the variation of the 
result with respect to the spectral index (0.5 x 10~ 9 ), and 0.6 x 10~ 9 from the uncertainty 
in the detector corrections (see Boone 2002, for details). Note that this result represents 
an average for the 2001 flare. Though spectral variability was reported at higher energies 
(Krennrich et al. 2002; Aharonian et al. 2002), STACEE-48 was not sufficiently sensitive to 
resolve this variation. 

Figure 3 is a plot of the correlation between the STACEE integral flux above 140 GeV 
and corresponding data taken by the Whipple 7-ray telescope (Mt. Hopkins, AZ, USA) 
above 390 GeV (Holder et al. 2001). Points represent data from both experiments whose 
time stamps agree to within 0.02 days (~29 minutes). A linear fit to the data yields a flux 
ratio of 4.3 ± 1.4. Assuming a simple power law spectrum for Markarian 421, this ratio 
corresponds to a differential spectral index of a — 2. 4+^4. If we include an exponential 
roll-off at E a = 4.3 TeV, the derived index drops to a — 2.3+®'^ both of which are consistent 
with the value of a = 2.14 from Krennrich et al. (2001). 

Assuming a differential spectrum of index a = 2.1 over STACEE's energy range, and 
after correcting for variations in the detector configuration, as well as for the observed star 
effect, we find the integral flux above 140 GeV to be (8.0 ±0.7 stat ±1.5 sys ) x 10~ 10 cm _2 s _1 . 
The first error is the statistical error in the Markarian 421 observations, and the second 
reflects the systematic errors discussed above, of which the error in the star correction is the 
dominant contribution (±1.3 x 10~ 10 ). 

5. Temporal Analysis 

While there is a substantial excess from Markarian 421 in the entire data set, individual 
on-off pairs generally do not exhibit highly significant detections. By combining the obser- 
vations of a single night, we can improve the flux significance at the cost of some temporal 
resolution. To obtain the STACEE-48 light curve for Markarian 421, we averaged the rate 
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over all the observations of each night, and assigned a date which was a weighted average of 
the Modified Julian Date in the middle of each run. 

Figure 4 includes the STACEE-48 light-curve for Markarian 421 between early March 
and late May. Filled circles are the counts per minute, and error bars are statistical only. 
For comparison, contemporaneous RXTE quick-look results (provided by the RXTE All 
Sky Monitor team) are plotted above, and Whipple observations (Holder et al. 2001) are 
plotted below. All the experiments observe an increased flux in the latter part of March, 
falling to lower values in April. Note that the coincident gaps in both the STACEE and 
Whipple coverage correspond to dates near the full Moon, during which Cherenkov telescopes 
generally cannot operate. 

The STACEE fluxes appear correlated with both the RXTE and Whipple observations. 
In fact, a discrete correlation analysis (Edelson & Krolik 1988) is consistent with the STACEE 
and Whipple observations being correlated on day timescales. However, the limited statistics 
in the nightly STACEE observations currently make analysis on shorter timescales infeasible. 

6. Conclusions 

STACEE has detected the BL Lac object, Markarian 421, with high significance in the 
140 GeV energy band, a hitherto unexplored region of its spectral energy distribution. The 
STACEE observed flux is consistent with other VHE 7-ray observations, and the temporal 
evolution of the STACEE observations appears similar to observations in both the TeV and 
X-ray bands. 

Figure 2 indicates that the STACEE measurements occupy an important region in 
Markarian 421 's spectral energy distribution near the peak of the IC hump. Future mea- 
surements of the spectrum by STACEE, along with simultaneous data at X-ray and TeV 
energies should help to further constrain synchrotron and IC blazar models for this source. 
STACEE is currently operating with 64 channels, each equipped with a flash ADC, and 
should be able to obtain more detailed spectral information in the near future. 

We are grateful to the staff of the National Solar Thermal Test Facility for their excellent 
support. Thanks to Gora Mohanty, Jeff Zweerink, Tumay Turner, Marta Lewandowska, 
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the National Science Foundation (under Grant Numbers PHY-9983836, PHY-0070927, and 
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Fig. 1. — STACEE-48 detector response to several different power-law spectra. Spectra are 
normalized to the same incident integral flux above 50 GeV for the given differential spectral 
index, a. 
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Fig. 2. — High-energy spectrum for Markarian 421. The STACEE butterfly assumes differ- 
ential spectral indices between a = 2.00 and 2.20. Note that, though the STACEE sensitivity 
peaks at 140 GeV, the vertex of the butterfly is much higher. STACEE data are averaged 
over the flaring period, while the Whipple curves represent high (I), medium (V), and low 
(VIII) flux levels over the course of the flare (see Krennrich et al. 2002, for details). EGRET 
data are from cycles 1 through 5, and thus represent an average low state. 
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Fig. 3. — Correlations between contemporaneous STACEE and Whipple observations from 
March to May of 2001. Points represent STACEE and Whipple runs whose time stamps 
match to within 0.02 days (29 minutes). The y 2 /doi for the fit is 26.1/22. The slope of 4.27 
is consistent with a differential spectral index of a = 2.14, given the STACEE and Whipple 
thresholds. 
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Fig. 4. — Markarian 421 light curve for Spring, 2001. RXTE observations are from the 2-10 
keV band, STACEE data are from the 50-300 GeV band, and Whipple observations span 
the 0.25-8 TeV band. References cited in the text. 



